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Abstract. The influence of extended convective mixing (over- 
shoot) on asymptotic giant branch stellar evolution is in- 
vestigated in detail. The extended mixing is treated time- 
dependently, and the efficiency declines exponentially with the 
geometric distance from the convective boundary. It has been 
considered at all convective boundaries, including the He-flash 
convection zone in the intershell region which forms during the 
thermal pulses. 

Both the structural and the chemical evolution are affected 
by the inclusion of overshoot. The main results include a very 
efficient third dredge-up which leads to the formation of car- 
bon stars of low mass and luminosity. A ^^C pocket which 
may serve as a neutron source for the s-process can form af- 
ter the third dredge-up has reached into the ^^C rich intershell. 
Overshoot applied to the pulse-driven convective zone during 
the He-flash leads to a deeper penetration of the bottom of this 
convective zone into the C/O core below the He-burning shell. 
This in turn causes *Hc to be less abundant in the intershell 
while ^^C and ^^O are more abundant compared to calcula- 
tions without overshoot. We show that overshoot at the He- 
flash convection zone as well as at the base of the envelope 
convection enhance the efficiency of the third dredge-up. Char- 
acteristic properties of the structural and chemical evolution of 
AGB stars are presented. 
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1. Introduction 

Towards the end of their lifetime, stars of low and intermediate 
mass (M < 8 Mq) evolve along the Asymptotic Giant Branch 
(AGB) stage ( |Iben and Renzini, 198"3| ; [Habing, 199^ ; Lattanzio 
and Boothroyd, 1997). The core of carbon and oxygen is sur- 
rounded by a sandwich like structure consisting of a helium 
burning shell, the hydrogen burning shell and the intershell re- 
gion in between. Evolved AGB stars undergo recurrent thermal 
instabilities of the helium burning shell (He-flash) which are 
refer red to as thermal pulses (TP) ( Schwarzschild and Harm, 
1965; [Weigert, 1966| ). Locally, helium burning peak luminosi- 
ties of logLne/ L© — 5 ... 8 are released and cause complex 



convective mixing events. The He-flash causes a pulse-driven 
convective zone (PDCZ) in the intershell region. After the He- 
flash the bottom boundary of the convective envelope may en- 
gulf the deeper regions where material previously synthesized 
by hydrogen and helium burning is present (third dredge-up, 
TDUP). 



Despite many s tudies of AGB evolution ( [ben, 1976 
[Schonberner, 1979|; ^^attanzio, 1986; Booth r oyd and Sack 



mann, 1988a; [Vassiliadis and Wood, 1993| ; [Blocker, 1995 



D'Antona and Mazzitelli, 1996; Forestini and Charbonnel 



1997; Straniero et al., 1997; Wagenhuber and Groenewegen, 



1998), important details, like the surface enrichment with nu- 
clear processed material from the deep interior (dredge-up) or 
the origin of ^^C which is an important source of neutrons 
for the synthesis of heavy elements in AGB stars, are not very 
well understood. To improve this situation we present the struc- 
tural properties as well as the chemical evolution of the inte- 
rior of AGB stellar models with convective overshoot. Convec- 
tive motions of matter approach the convective boundary with a 
non-zero velocity and penetrate into the formally stable region. 
These overshooting flows lead to extra mixing of elements. 



Canuto (1998) has pointed out that overshoot is a dynamical 



consequence of Newton's laws and as such is unavoidable. Pre- 
vious studies have shown that models with overshoot can ac- 
count for several observed properties of AGB and post-AGB 
stars (Blocker et al., 1997; Herwig et al., 1997, 1998, 1999b). 

We give a brief review of AGB star properties relevant for 
our new models (Sect. some remarks on the stellar evolution 
code as well as some information on overshoot and its treat- 
ment in our models (Sect.||). We explain the differences be- 
tween models without overshoot and with overshoot by looking 
separately at the two relevant convective boundaries: the bot- 
tom of the envelope convection (Sect.^) and the bottom of the 
PDCZ (Sect.^). The chemical and structural surface properties 
of models with overshoot are described in Sect. ^. Conclusions 
are presented in Sect. ^ 

2. Some properties of AGB stars, stellar models, and 
overshoot 

Carbon stars and the third dredge-up AGB stars show a large 
range of surface C/O ratios from about the solar value ~ 0.4 
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up to well above unity. In particular many low-mass AGB stars 
show up as carbon stars with C/0 > 1 (Smith et al., 1987; 
Frogel et al., 1990). Half of all planetary nebulae are carbon 



rich as well (Zuckerman and Aller, 1986). Already the early 
stellar evolution calculations ([ben, 1975; Iben, 1977; Iben and 



Truran, 1978; [Sackmann, 1980| ) pointed towards the TDUP to 
provide a link between the intershell region where carbon is 
present and the bottom of the envelope convection. ^ How- 
ever, these models have not shown sufficient dredge-up for 



AGB models with low (core) masses. Lattanzio (1989) used 
a method to determine the convective boundary which ensures 
that the ratio of radiative and adiabatic temperature gradient 
approaches unity smoothly at the convective boundary. This 
method favors the occurrence of the TDUP even for low mass 



stars of solar metallicities. Hollowell & Iben ( |1988| ) found 
for low mass and low metallicity that the dredge-up of pro- 
cessed material occurs if some additional mixing is assumed. 
Wood ( |1981| ) and Boothroyd & Sackmann ( |l988b[ ) showed that 
lower metallicity and increased mixing length parameter en- 
hances TDUP. Straniero et al. ( |1997| ) find the TDUP for low- 
mass AGB stars and relate its occurrence to increased numer- 
ical resolution. Calculations which do not find dredge-up for 
low core masses are the rule rather than the exception (Vassil- 
iadis and Wood, 1993; Blocker, 1995; Forestini and Charbon- 



nel, 1997; Wagenhuber and Groenewegen, 1998; Langer et al., 
1999). Herwig et al. ( |1997[ ) presented first results of a 3 M© 
stellar model sequence which had been computed with con- 
vective overshoot, treated time-dependently at all convective 
boundaries. With this approach TDUP was very efficient for 
low core masses and solar metallicity. 

Intershell abundances The abundance distribution in the in- 
tershell is of great importance because in this region the major 
nuclear burning and mixing processes associated with the He- 
flash and the TDUP take place. The intershell abundances are 
determined by nuclear burning during the quiescent interpulse 
phase (mainly hydrogen burning), the nuclear processing dur- 
ing the TP, the convective mixing in the pulse-driven convective 
zone, and third dredge-up. The intershell abundance is signif- 
icantly affected if overshoot is applied to the PDCZ because 
additional material is mixed into the intershell from the C/0 



core (Sect. 5.2) 



It affects the local nuclear production, the surface enrich- 
ment and also the structure of the star. It is also of relevance for 
the interpretation of surface abundances of hydrogen-deficient 
post-AGB stars (see e.g. Schonberner, 1996; Iben et al., 1996 
or Werner etal., 1999). 



Nei ther Schonberner ( 1979 ) nor Boothroyd & Sackmann 
( 1988a) considered overshoot and both found that the abun- 



dance distribution in the intershell evolves with each TP. Start- 



The third dredge-up is to be distinguished from the first dredge-up 
which occurs during the first ascent of the giant branch after central 
hydrogen exhaustion and the second dredge-up which occurs during 
the early AGB phase for more massive AGB stars after the end of core 
helium burning. 



ing with an almost pure He composition from previous hydro- 
gen burning at the first TP, the abundances approach typical val- 
ues of (He/C/O)=(0.76/0.22/0.02) (mass fractions) after a few 
TPs. 

^^C production The observed correlation of s-process ele- 
ments and carbon in low-mass stars ( Smith and Lambert, 1990| ) 
points to low-mass stars as a likely site for n-capture nucleosyn- 
thesis ([GalUno et al., 1997| ; |Wallerstein et al., 199*^ . GalUno 
et al. (1998) have demonstrated that, hydrogen ingestion 



from the envelope into the intershell region at the end of the 
TDUP is assumed then the subsequent formation of ^^C via the 
nuclear reaction ^^C(p, 7)^^N(e"'"i/)^^C does indeed predict a 
neutron exposure which leads to a s-process enhancement in 
compliance with the solar main component of heavy elements. 
However, the physical mechanism of H-ingestion is unclear. 
Depth-dependent overshoot naturally predicts the formation of 
a ^^C pocket as needed for the s-process nucleosynthesis. 



Mixing and overshoot Mixing of elements in stars is at- 
tributed to a number of processes ( ^insonneault, 1997 ) of 
which convection is the most effective. Stellar evolution mod- 
els commonly employ the mixing-length theory (MLT) (Bohm- 
Vitense, 1958) or some descendent thereof. The boundary of 
convective instability is determined by the local Schwarzschild 
condition Viad > Vad where Vad and Vrad are the adia- 



batic and radiative gradient (Kippenhahn and Weigert, 1990). 



However, neighboring layers are related by inertia, momentum 
transfer and the equation of continuity and therefore convective 
elements might overshoot beyond the boundary of convection 
(e.g. Shaviv & Salpeter, 1973; Maeder, 1975; Roxburgh, 1978; 
Bressan et al.,1981; Langer, 1986). No commonly accepted 
quantitative the oretical descrip tion of convective overshoot 
currently exists ( l^enzini, 1987 ). 



Comparison of stellar models with observational findings 
indicates that convective overshoot occurs in real stars (Mer- 



milliod and Maede 
dersen et al., 1990 


r, 1986; 


Maeder and Meynet, 1989 


; An- 


IStothers, 1991; 


Napiwotzki et al., 1991 


A.longi et al., 1991; 


Alongi et al., 1993 


; [Schroder et al., 1997 


fCozhurina-Platais et al., 1997). Overshooting in stellar evolu- 



tion calculations has been simulated by an extension of the in- 
stantaneous mixing beyond the convective boundary (instanta- 
neous overshoot). The widely used set of models by Schaller 
et al. (1992) has been calculated with an overshoot distance of 
one fifth pressure scale height (aover — 0.2) which was cali- 
brated in order to fit the observed terminal age main sequence 
of 65 stellar clusters and associations. Improving on the instan- 
taneous treatment, Deng et al. (1996a; 1996b) and Salasnich 



et al. (1999)) have explored the effects of turbulent diffusion in 
the overshoot region of massive stars. 

Due to the insights of two- as well as three-dimensional 
hydrodynamical simulations convection is nowadays pictured 
in terms o f downdrafts and up-flows rather than as a hierarchy 
of eddies ( Stein and Nordlund, 1998 ). Two-dimensional simu- 
lations showed that prominent downward-directed plumes can 
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overshoot a substantial distance into the stable region (Hurl- 
burt et al., 1986; ^urlburt et al., 1994| ; [Freytag et al., 1996| ). 
These models show emerging and vanishing patterns of curls, 
fast narrow downdrafts and broad up-flow regions. The turbu- 
lent velocity field decays exponentially beyond the convec- 
tive boundary (see also Xiong, 1985; Asida & Arnett, in prep.). 
These results have been applied via a time-dependent treat- 
ment of convective mixing to low- and intermediate mass 



stellar models by Herwig et al. (1997; 2000) (low mass TP- 
AGB), Ventura et al. (1998) (main sequence) and Mazzitelli 



et al. ( I1999D and Blocker et al. (gOO0|) (massive AGB stai's, 
HBB). 



3. The stellar evolution code and computations 

General aspects of the code The results presented in this pa- 
per have been obtai ned w ith the stellar evolution code de- 
scribed by Bl ocker ( 1995 ) with the modifications given in 
Herwig et al. ( 1997 ). Here, we will add only a few relevant 
details. The C (a, 7 V^O reaction rate has been taken from 
Caughlan & Fowler ( 1988 ) and multiplie d by a factor of 1.7 
as recommended by Weaver & Woosley (1993). The reaction 



rates ^^0(p,7)^^F and ^^0(p, a]^^ have been taken from 



the co mpilation of Landre et al. ( 199C ) as described in El Bid 
( 1994 ). The model sequence without overshoot has been com- 
puted with the latest OPAL opacity tables (Iglesias and Rogers, 
1996). This choice had purely technical reasons but does not 
affect the comparison. 

The adjustment of the numerical resolution (geometry and 
time) plays an important role in the computation of AGB mod- 
^rost and Lattanzio, 1996 ; 



els (Straniero et al., 1997 



Mowlavi, 

1999). The resolution is adjusted from model to model ac- 
cording to the changes of the state variables as a function 
of time and location within the star For example, the lumi- 
nosity due to nuclear processing of helium Lhc should not 
change more than 5% between two models. The models have 
~ 1900 . . . 2300 mass grid points which accumulate around the 
core-envelope interface where the pressure and density drop 
steeply. The model time steps range from about a dozen years 
during the quiescent H-burning interpulse phase down to typ- 
ically a few days during the He-flash and the following phase 
where dredge-up may occur (Sect. 4.1). 



Exponential diffusive overshoot For the time-dependent treat- 
ment of ove rshoot mixing we follow the prescription of Freytag 
et al. ( 1996 ). The particle spreading in the overshoot region can 
be described as a diffusion process and the diffusion coefficient 
Dqy can be fitted by a formula like 



(1) 



where z denotes the geometric distance from the edge of the 
convective zone and <c is a characteristic time scale for the con- 
sidered convection zone. Combining this with the finding of an 



exponential decay of the velocity field, Freytag et al. ( 199^ 
give the diffusion coefficient in the overshoot region as 



DoY = Do cxp 



-2z 

~h7 



P' 



(2) 



where is the velocity scale height of the overshoot con- 
vective elements at the convective boundary. can be ex- 
pressed as a fraction / of the pressure scale height Hp. The 
quantity Dq is the expression fedge^c in Freytag et al. (1996: 
Eq. 9). Here we have approximated Dq by the diffusion coef- 
ficient in the convectively unstable region near the convective 
boundary redgc- Note that Dq is well defined because Dc ap- 
proaches the convective boundary with a small slope and drops 
almost discontinuously at rodgo- The diffusion coefficient Dc in 
the convection zone can be derived from the MLT (see Langer 
et al., 1985). In the radiative zone where the diffusion coeffi- 
cient has dropped below -Dqv'' — 1 ' lO^^cm^/s no element 
mixing is allowed (Z^iad = 0). 

The free parameter / in Eq. (||) describes the efficiency of 
the extra diffusive mixing. For larger / the extra partial mixing 
beyond the convective edge extends further. In their simulations 
Freytag et al. ( 1996 ) found velocity scale heights of the order 
of a pressure scale height (/ ~ 1). However, under adiabatic 
conditions the convective flows move so fast that they will not 
be able to exchange heat efficiently with the surrounding when 
entering into the stable region. Accordingly a stronger deceler- 
ation of the plumes is likely in the stellar interior compared to 
the situation of the shallow surface convection zones. In order 
to obtain an order of magnitude for / in the stellar interior it 
has been scaled to reproduce the observed width of the main 
sequence. A parameter of / ~ 0.016 was found to reproduce 
the models of Schaller et al.(|l992|) (see end of Sect.0). Us- 



ing AGB overshoot models as a starting point for post-AGB 
models of H-deficient objects has confirmed this assumption 
a posteriori, at least for the bottom of the He-flash convection 
zone ( ^erwig et al., 1999a[ ). 

The overshoot prescription (Eq. (||)) has been applied irre- 
spective of /i-gradients which may decrease the overshoot ef- 
ficiency (Langer et al., 1983). The analysis of old open clus- 
ters even suggests that overshoot may not expand against /i- 
gradients in low mass stars (Aparicio et al., 1991). Theoreti- 
cally, Canuto (1998) has recently argued that /^-gradients do 
decrease the overshoot efficiency, but do not prevent the effect 
entirely. In our study, this aspect becomes important at the bot- 
tom of the envelo pe c onvection during the third dredge-up. As 
explained in Sect. the results are not significantly dependent 
on the amount of overshoot at this convective boundary. 

The stellar evolution computations The findings in this work 
are based on two model sequences of 3 Mq and 4 Mq with 
overshoot. For comparison, one 3 Mq TP- AGB sequence has 
been computed without overshoot. It has been started from a 
3 Mq overshoot model taken just before the occurrence of the 
first TP (Fig.^ ^, |o[ |l5]). In order to study the dependence 
of the evolution of the intershell abundance as a function of the 



4 



F. Herwig: AGB evolution with overslioot 




(a) TPS, no OV 



PDCZ 



h- . / 1 I- I . I . I . I ■ I ■ I . I . I .H 



-1 1 



2 3 4 5 6 
t/[100 yr] 



op 
O 



CL 
H 



I ' I ' I ' I 



I ' I ' I ' I ' I 



'CB 
Mu 



(b) TP14, no OV 



h- /. 1 1 ■ I . I . I . I . I ■ I . I . I . H 



-1 01234567 
t/[100 yr] 



I ' I 



-1 1 



I ' I ' I ' I ' I ' I ' I 

McB 

Mh 



(c) TP14, DUP: OV 



h- . I ■ I . I . I . I ■ I ■ I . I . I 



3 4 5 
t/[100 yr] 



CO 

o 



CL 
H 



I ' I ' I ' I ' I ' I 

McB 

Ml 



(d) TP14^DUP:0V 
PDC^: OV 




3 4 5 
t/[100 yr] 



Fig. 1. Evolution of mass coordinate of convective boundary (Mcb, solid line) and H-free core (Mh, dashed line). The mass 
scale has been set to zero at the mass coordinate of the H-free core at the time of the He-flash (TPS: Mu = 0.6421 M©; TP14: 
Mh = 0.6917 Mq). The time scale has bee set to zero at the peak of the He-flash luminosity (TPS; to = 227034 yr; TP14: 
to — 677076 yr). The spike (solid line) at t — Oyr shown in (a), (b) and (d) is the top of the pulse-driven convection zone of 
the He-flash. The upper solid line shows the bottom boundary of the envelope convection zone. Panel (a): 5**^ TP of the 3 Mq 
sequence without overshoot; panel (b): 14*'' TP of the same sequence; panel (c): the same TP as in (b) started after the PDCZ 
has vanished and with overshoot at the bottom of the envelope convection zone; panel (d): again the 14*'' TP recalculated from 
before the onset of the He-flash, with overshoot at the PDCZ and during the following TDUP episode. 



parameter / a few more sequences with different /-values have 
been computed from this starting model (Fig. |l2|) and followed 
over three thermal pulses. In addition, the eighth TP cycle of 
the 3 Mq overshoot case has been recomputed for different / 
values (Fig. [l^, |l^. In Fig. |l] a recomputation of the 14"^ TP of 
the 3 Mq sequence initially without overshoot is shown. 

A solar-Uke initial composition, (Y, Z) — (0.28, 0.02), has 
been chosen corresponding to Anders & Grevesse (1989). The 
mixing-length parameter of the MLT is q;mlt = 1- 7 which has 
been calibrated by reproduction of the solar surface parame- 
ters. A Reimers-type mass loss with efficiency rj = 1 has been 
applied. The Reimers mass loss formula has not been designed 
for AGB stars, however for the purpose of this study mass loss 
is not important. The 3 Mq (4 Mq) model sequence has a to- 



tal mass of 2.686 Mq (3.689 Mq) at the 12*'' (11**') TP (see 
Tabled 

Overshoot has been applied during all evolutionary stages 
and to all convective boundaries. Overshoot applied to the main 
sequence core convection has indirect impUcations for AGB 
star models as it shifts the core masses at the first thermal 
pulse upward and decreases the progenitor mass (A/up) limit 
for C/O white dwarfs (/Vlongi et al., 1993). If overshoot is also 
applied during the second DUP phase this effect is counterbal- 
anced somewhat and M^p — 6.5 (T. Driebe: priv. com., see 



also Blocker et al. (200C) and Weidemann (2000)). This is just 



above the value obtained from models in which overshoot is 
only applied during the main sequence phase (B ressan et al., 
1993). 
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Fig. 2. Temperature gradients, opacity and abundances for the 
mass fraction in the AGB star during the TDUP (/ — 0.016) 
after the eighth TP with overshoot, t — 216 yr after peak he- 
Uum luminosity. Compare Fig. ^ there the bottom of the con- 
vective envelope does not reach into the He-rich buffer zone at 
Mr = 0.6529 Mq, whereas in the case shown here the enve- 
lope has started to engulf the He-rich buffer zone. 



4. Overshoot at the bottom of the convective envelope 

Overshoot at the bottom of the convective envelope has an ef- 
fects on the modeling of the TDUP and on the formation of a 
^•^C pocket after the end of the dredge-up phase. 



4.1. The third dredge-up 

It has been a matter of discussion whether convective overshoot 
at the botto m of the conve ctive envelope is a prerequisite for 
dredge-up ( ^owlavi, 1999 ). This question is closely linked to 
the general problem of determining the convective boundary 
of mixing within the framework of a local theory of convec- 
tion together with its numerical implementation. As reported 
by previous authors (e.g. Paczyfiski, 1977), our models without 
overshoot develop a discontinuity in Viad when the bottom of 
the envelope convection zone reaches the He-rich buffer layer 
This discontinuity can prohibit the TDUP under certain condi- 
tions. In our models of AGB stars with not so massive cores the 
flux on the stable side of the convective boundary is not large 
enough to lift Vrad above Vad ■ 

In Fig. |l] the situation in the dredge-up region is shown for 
several thermal pulses computed with different assumptions on 
overshoot. After an early TP of the 3 Mq sequence computed 
without any overshoot (panel a) the bottom of the envelope 
convection does not approach the H-free core close enough for 
any dredge-up. For more advanced TP (panel b) the convection 
clearly approaches the H-free core for a period of about 300 yr 
but due to the effect of the abundance discontinuity shown in 
Fig. |] the lack of overshoot prohibits the TDUP. This is evi- 



dent from panel (c) which shows the same situation with over- 
shoot switched on below the bottom of the envelope convec- 
tion. In this case a considerable dredge-up of material is pos- 
sible. Panel (d) shows that overshoot applied also to the pre- 
ceding PDCZ increases the DUP efficiency even further, which 
is discussed in detail in Sect. ||. Note that if overshoot only at 
the base of the convective envelope had been applied for the 
^th gjjQ.^jj jjj Fig.|l], this would have had no effect on DUP 
because the stable layer between the bottom boundary of con- 
vection and the H-free core is too large compared to the tiny 
extent of the overshoot layer. Fig. ^ demonstrates the situation 
of panel (b) of Fig. |l] as a set of profiles. The bottom of the con- 
vective envelope recoils at the He-rich region which encloses 
the core, despite the fact that contact between the convective 
envelope and the H-free core is established over a few hundred 
years. In Fig.|^ the time after the eighth TP of the 3 Mq se- 
quence with overshoot has been recalculated without overshoot 
at the bottom of the envelope convection. A profile at the begin- 
ning of the TDUP of the corresponding original sequence with 
overshoot is shown in Fig.^. Close inspection and compari- 
son with Fig.|| shows that 10^^ Mq of He-rich material have 
already been dredged-up. The additional mixing removes the 
sharp abundance discontinuity. This in turn dissolves the sharp 
discontinuity in the radiative gradient Vrad which sensitively 
depends on the opacity. In the envelope, hydrogen has a mass 
fraction of ~ 0.70 and a higher opacity than helium. If this ma- 
terial is mixed into the stable region the radiative gradient will 
be lifted and dredge-up can occur 

This principle is not dependent on the efficiency of over- 
shoot within a factor of maybe two (Fig. |l3|). At the lower edge 
of the AGB envelope the convective velocities are of the or- 
der of ~ Ikm/s while very efficient mixing at this convec- 
tive boundary is still achieved by convective velocities down to 
< Icm/s. An instantaneously mixed region can be identified 
in which the downward directed flows decelerate from the ini- 
tial velocity at the edge to the much lower velocity where the 
chemical abundance is not homogenized. Here, the turbulent 
velocity is so slow that mixing becomes inefficient. If the over- 
shoot efficiency is small the region of homogenized composi- 
tion may be very thin. This will still suffice to allow TDUP if 
the resolution is chosen such that overshoot is numerically able 
to smear out the abundance discontinuity (see below). Fig.^ 
gives an impression of the masswise size of the overshoot re- 
gion which can be barely identified in the abundance profiles 
on this scale. If the time steps are too large with a small over- 
shoot, the smoothing effect on the abundance discontinuity is 
not efficient and DUP is prevented as is the case in the exam- 
ples shown in Fig. [l|b and Fig. ^. 

In our calculations without overshoot Vrad and Vad are 
first compared at each mass grid zone which defines stable and 
unstable grids. If a discontinuity has already built up, grids are 
included at the discontinuity in order to still resolve the abun- 
dance gradient. No mixing is allowed from an unstable grid to 
a stable grid. This approach closely realizes the MLT which 
predicts that no convective bubble can leave the convectively 
unstable region. However, one can argue equally well that the 
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Fig. 3. A time series of the region of the former hydrogen-burning shell after the eighth TP of the 3 M© sequence. The dredge- 
up phase after the thermal pulse of the overshoot sequence has been recalculated without overshoot and as a result dredge-up 
is prevented (compare Fig.|l]) and Fig.|[ Abundances in mass fractions of "^He (short-dashed), ^^C (dot-dashed) and hydrogen 
(dashed) refer to the left scale while the adiabatic (Vad, long-dashed) and radiative (Vjad, solid line) temperature gradient refer 
to the right scale. For each panel the abscissa ranges from the bottom of the envelope at the right to the top of the intershell region 
at the left. In the middle of each panel a region of almost pure ''He has been formed by hydrogen-burning during the previous 
interpulse phase. The lower boundary of the He-rich region is formed by the He-flash convection zone, which has at < = Oyr 
homogenized the intershell layers below = 0.6527 Mq. The four snapshots show how the bottom of the envelope convection 
zone advances downward by mass during the phase when TDUP would occur with overshoot. The order is from top to bottom 
and left to right and the ages are given with respect to the flash peak of the previous pulse. Note how the discontinuity of Vrad 
is related to the abundance discontinuity which starts to build up (t — 252 yr). The discontinuity grows further (t — 398 yr) 
and finally when the region resumes contraction the bottom of the envelope convection retreats without being able to dredge-up 
material from the H-exhausted regions. 



average of the gradients between neighboring grids should be 
compared. This criterion would allow material in stable grids 
to be mixed with material from the neighboring grid if the lat- 
ter is unstable. Under certain conditions dredge-up may occur 
with the second treatment which would not occur with the first 
treatment. 

During the dredge-up the upper part of the intershell region 
is engulfed by the convective envelope at the rate A/dup = 
dnicdgc/dt where modgc is the mass coordinate of the con- 
vective boundary. Then, the extent AA/h (in Lagrangian co- 
ordinates) of the stable layer which is practically homogenized 
by extra mixing, depends on the overshoot efficiency and the 
rate Afoup- AAfh is larger for a larger efficiency. Moreover, 
AAfh will be larger if A/dup is smaller and vice versa be- 
cause element mixing is treated time dependently. With the ef- 



ficiency / = 0.016 and A/dup - 2 • lO^^Mo/yr we find 
AA/h — 10^^ Mq . In order to resolve the overshoot region the 
mesh size at and around the convective boundary should not 
exceed AA/h/10. 

The time step is limited by the extent of the homogenized 
stable layer Current stellar evolution calculations typically sep- 
arate the solution of the structural equations and the solution of 
the mixing equations (often a diffusion-like equation) in order 
to save computational resources. This separated solution is jus- 
tified if the two processes described by the two sets of equa- 
tions operate on different time scales. If the numerical time 
step is of the order At > AA/h/ A/dup then the mass co- 
ordinate of the envelope convection bottom mcdgc after this 
time step would be below the previously homogenized layer: 
m-cdgc(i + At) < rricdgcit) — AA/h. In that case the sepa- 
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rated solution of the structure equation and the equations of 
abundance change is no longer valid because the abundances 
Y{t + At) at mcdgc(^ + A<) cannot be approximated by the 
abundances Y(t). Only if this approximation is valid can the 
separate computation of structure and abundance change give 
correct results. This is only the case if At < AA/h/A/DUP- 

The actual spatial resolution is controlled by a combina- 
tion of criteria. Throughout the star we insert mass shells if 
any of the conditions AlogL > 0.01, AlogP > 0.02, 
AlogX(i2C) > 0.02 or AlogX(4He) > 0.02 is fulfilled 
between two mass shells. After the He-flash, when the dredge- 
up is expected we increase the resolution in a region of about 
0.06 • Mcore around the mass coordinate at the bottom of the 
envelope convection zone by requiring AlogP < 0.01. This 
prescription for the spatial resolution ensures the above men- 
tioned mesh size constraints in an adaptive way. Any further 
increase of resolution does not affect the amount of dredge-up 
obtained. Adapting the mesh may cause additional numerical 
diffusion. In the situation of the advancing convective bound- 
ary associated with the TDUP this does not prevent the discon- 
tinuity (Fig.^ if no minimum grid size is enforced. However, 
numerical diffusion may easily distort the abundance profiles 
in the ^^C pocket (Fig. ^ when the grid density after the end of 
the TDUP is first reduced and then increases afterwards when 
H-burning starts. 

4.2. The formation of the ^^C pocket and of an adjacent ^^N 
pocket 

During the TDUP the hydrogen rich and convectively unsta- 
ble envelope has contact with the radiative carbon rich in- 
tershell region. If the boundary of the envelope convection 
zone is treated according to the exponential diffusive overshoot 
method hydrogen diffuses into the carbon rich layer below. 
When the bottom of the envelope convection zone has reached 
the deepest position by mass, a region forms were protons from 
the envelope and ^^C from the intershell region coexist (Fig.^ 
top panel). Note that the whole mass range shown in Fig. ^ is 
part of the overshoot region of the convective envelope. The 
solid line in the top panel shows the exponentially declining 
diffusion coefficient at the end of the dredge-up episode. In the 
right half of the panel D is large enough to cause a homoge- 
neous element mixture. The profiles have developed cumula- 
tively over a period of w 50 yr during which the mass coor- 
dinate of the base of the envelope convection remains almost 
constant. In the left part the hydrogen abundance decreases as 
^^C increases. Thus, the H/^^C number ratio decreases con- 
tinuously from ^ 750 in the envelope to zero in the intershell 
region. 

When the temperature increases during the succeeding evo- 
lution proton captures transform ^^C and H according to the 
continuously decreasing H/^^C ratio into a ^^C pocket and 
a prominent ^"^N pocket side by side. The ^^C forms in the 
deeper layers where the H-abundance is below ~ 5%. Here 
the newly formed ^"^C cannot be further processed into ^''N 
because all protons are akeady consumed. The ^^C abun- 



dance reaches a maximum mass fraction of 0.09 at Afr = 
0.649144 M0. Above this mass coordinate the H/^^C ratio is 
larger and after the formation of some ^^C there are still pro- 
tons available which transform ^"^C into ^^N. The ^■^C/^^N ra- 
tio changes continuously in this region and at any position it is 
almost entirely determined by the initially present H/^^C ratio. 
Moreover, a ^"^N pocket is a direct consequence of the contin- 
uous variation of the H and ^^C abundances from the core to 
the envelope. Therefore, if the efficiency of mixing is decreas- 
ing smoothly in the overshoot region a substantial ^''^N pocket 
will inevitably form. The maximum abundance is given by the 
amount of ^^C in the intershell (see Sect. 5^). In a narrow re- 
gion (albeit larger than the region occupied by the ^"^C pocket) 
^^N becomes the most abundant isotope at a mass fraction of 
0.45. The bend of the profile, e.g. at A/,. = 0.649146 M© orig- 
inates from the changing abundance ratios of ^^C, H and ^"^N 
which are relevant for the production of ^^C. 

Note that at the time when proton captures start to form 
the ^■^C and ^"'N pockets the radiative region above the burn- 
ing region is well established. The isotopes made in this region 
cannot reach the surface at this time. Instead they are processed 
under radiative conditions ( Straniero et al., 1995 ) and the prod- 
ucts of this processing will be engulfed by the next He-flash 
convection zone. Thus, the production of ^■^C after dredge-up 
with overshoot does not decrease the observed ^^C/^^C ratio in 



giant stars as suspected by Wallerstein & Knapp ( |1998[ ). 

During the whole interpulse period the region displayed in 
Fig.^ is heating and contracting (Fig.^). When the tempera- 
ture has reached about T = 10^ K, is destroyed (Fig.^ 
by ^^C{a,n)^^0 and neutrons are released. Towards the end 
of the interpulse period (third panel Fig.^ most of the ^^C is 
burnt by a capture. The remaining small fraction is engulfed 
by the outwards reaching He-flash convection zone of the ninth 
pulse. This is displayed in the bottom panel of Fig.^. The 
straight solid line shows the diffusion coefficient in the over- 
shoot region at the top of this convection zone. Although the 
absolute value of D is larger than in the top panel, the impact 
on the abundance appears to be smaller. This effect is caused 
by the smaller velocity of the convective boundary of the enve- 
lope during the dredge-up compared to that of the boundary of 
the He-flash convection zone. 

Fig.^ shows that after more advanced pulses (like the 
eighth in this case), the ^'^C formed at the hydrogen - carbon 
interface d oes indeed burn almos t completely under radiative 
conditions ( [Straniero et al., 1995 ). The ^■^C which formed af- 
ter earlier pulses (like the fifth in this case) burns only partly 
under radiative conditions. After early thermal pulses a certain 
fraction of the ^^C will be mixed into the pulse-driven con- 
vective zone of the succeeding TP where it is processed under 
convective conditions in the intershell. In Fig. ^ the density and 
temperature at the position of the ^'^C pocket are given for the 
case of the eighth TP. This shows the conditions under which 
the neutron capture operates. 

With our overshoot description and assumption about the 
efficiency (/ = 0.016) the layer in which the ^■^C pocket forms 
is very thin (see Fig.0). The total mass of ^"^C contained in the 
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Fig. 4. Development 
and destruction of the 
^^C pocket after the 
eighth TP of the 3 Mq 
sequence. The panels 
show the abundance pro- 
files of H/^C and other 
isotopes in a mass range 
that contains the contact 
region of envelope and 
core at the end of the 
TDUP. The top panel 
shows a model after the 
end of the dredge-up 
episode (about 500 yr 
after the He-flash). The 
second panel shows a 
model about 1800 yr 
later when hydrogen 
burning has set in again. 
The third panel shows 
a model at the very end 
of the interpulse phase 
when ^•^C has already 
been destroyed by re- 
action ^^C{a,n)^^0. 
The last panel shows the 
profiles at the onset of 
the next thermal pulse. 
The solid lines in the 
top and bottom panel are 
the diffusion coefficient 
(right scale) of these 
models. In the top panel 
the diffusion coefficient 
results from the over- 
shoot of the envelope 
convection while in the 
bottom panel it is the 
diffusion coefficient of 
the upper overshoot 
zone of the He-flash 
convection zone during 
the ninth pulse 



pocket comprises only « 2 ... 4 • 10"'' Mq. Estimates of the 
s-process element distribution from the ^■^C neutron source de- 
mand that the pocket contains w 1 . . . 2 ■ 10^^ in order t o 
model the main s-process component ( |Straniero et al., 1995 ). 
However, the numbers are not comparable at face value because 
the result of Straniero et al. is based on models without over- 
shoot. Their hydrogen profile in the carbon rich region is not 
the same as our overshoot hydrogen profile. Also, the contri- 



bution and modification of the s-process element distribution 
from the ^^Ne(Q;, n)'^^Mg neutron source reaction during the 
high-temperature phase of the flash should be different in our 
models (see Sect. 5A and^. Finally, the additional ^^N pocket 
which forms due to the overshoot model must be considered. 
On one hand ^''N is known to be an important neutron poi- 
son because of its large {n,p) cross section. On the other hand 
a large fraction of the neutrons lost by this reaction may be 
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Fig. 5. The maximum ^'^C abundance of the pocket after the 
fifth and eighth pulse of the 3 Mq sequence displayed over one 
pulse cycle respectively. The zero time is set to the flash peak 
respectively. 
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Fig. 6. The temperature and density at the location of the maxi- 
mum ^•^C abundance for the eighth TP corresponding to Fig.||. 



reproduced by ^^C(p, n) i^N(/3+ Y^C[a,nY^O. Moreover the 
additional amount of ^^N to be ingested into the next He-flash 
convection zone will also be converted into additional ^^Ne. 
All those processes will affect the s-process nucleosynthesis, 
which cannot be reduced to the mere amount of ^'^C produced. 

Preliminary tests with different overshoot efficiencies have 
shown that the shape of the abundance profiles in the region 
of the ^•^C pocket is scaled with respect to the mass coordi- 
nate but conserved with respect to the abundance ratios. Inter- 
vals with certain mass ratios of, e.g. ^■^C and hydrogen are just 
larger with a larger /. This means that a dedicated study of the 
s-process nucleosynthesis with a ^■^C pocket according to our 
overshoot description should find that within a certain range of 
efficiencies the f-value does only determine the total amount of 
s-process elements in the dredged-up material but not its distri- 
bution. It should be investigated whether the functional form of 
an exponential velocity decay can reproduce a s-process distri- 
bution in compliance with the solar main component. 

Apart from overshoot the process of rotationally induced 
mixing is a promising mech anism for the origin of the ^"^C 
pocket ( l^anger et al., 1999 ). It is surprising that the over- 
all amount of ^■^C found in the pocket due to rotation is al- 
most identical to the amount found here with overshoot (a few 

10-^ Mq). 



5. The role of overshoot at the boundaries of the He-flash 
convection zone 

During the He-flash the intershell region becomes convectively 
unstable due to the huge energy generation of the He-burning 
shell. Overshoot at the bottom of this He-flash convection zone 
leads to a deeper mixing from the region below the He-shell 
(the C/O core) into the intershell. This process could be called 



the fourth dredge-up ([ben, 1999) or intershell dredge-up. Both 



the abundance in the intershell and also the structural parame- 
ters in the convection zone are affected by overshoot at the base 
of the He-flash convection zone. In particular, the change of the 
intershell abundance and the larger en ergy generation durin g 
the TP leads to very efficient TDUP ( [Herwig et al., 1999b| ). 
If these overshoot AGB models are used as starting models 
for the post-AGB evolution the long-standing discrepancy be- 
tween observed surface abundances of H-deficient post-AGB 
stars and stellar models of this phase can be resolved (Herwig 
et al., 1999a). 

5.7. Intershell structure during the flash 

During the onset of the He-flash, when helium-burning in- 
creases rapidly, the layers just above the position of maximum 
energy generation and temperature become convectively unsta- 
ble and the pulse-driven convection zone develops (Fig.|^). In 
this situation overshoot at the bottom of the PDCZ supplies 
additional fresh helium into the underlying nuclear burning re- 
gion and supports the nuclear runaway. The situation is some- 
what similar to hot-bottom burning were the hydrogen burn- 
ing shell obtains additional nuclear fuel from the deep enve- 
lope convection and accordingly increases its energy produc- 
tion. Similarly the extra mixing below the PDCZ leads to a 
more violent He-flash with larger peak He-burning luminosi- 
ties (Fig.^ and larger mass AMpocz of the convectively un- 
stable region (Fig.|^d). The stronger He-flashes are responsi- 
ble for the structural and abundance consequences discussed 
below. The larger extent of the He-flash convection zone has 
immediate consequences, e.g. for the s-process nucleosynthe- 
sis: Heavy elements produced during the interpulse phase un- 
der radiative conditions are diluted in the convective intershell 
during the pulse before they are dredged-up into the envelope. 
For a larger He-flash convection zone the dilution effect will be 
larger. 

Apparently there is no noticeable effect of overshoot at the 
top of the PDCZ. For the metallicity and overshoot efficiency 
considered here the minimum mass layer remaining between 
the H-rich envelope and the top of the PDCZ at its largest exten- 
sion is independent on overshoot as can be seen from Fig. 

The sequence without overshoot does not show any TDUP, 
and accordingly the whole intershell with the He-burning, and 
the H-burning shell at each boundary respectively is shifted 
outward by mass. For the 3 M0 sequence with overshoot 
dredge-up sets in at the third TP and the effective core mass 
growth (AA/h — AA/dup) is slowed down (Fig. ^,b). At the 
last computed TP of the 3 Mq sequence the intershell no longer 
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evolves outward with mass. During this stationary shell burn- 
ing the nuclear fuel for both the hydrogen and helium shell is 
transported downward by dredge-up. Over the pulse cycle the 
dredged-down hydrogen is first converted into helium which 
is then further mixed down by the He-flash convection zone. 
Finally the burning products are exchanged by dredge-up with 
fresh envelope material. In this situation dredge-up is not only 
important to mix processed material up but also to mix fresh 
material down and feed the stationary shells. 

For the 4 Mq sequence the core mass actually decreases 
from pulse to pulse because A > 1 (Table ^. In that case the 
nuclear burning shells are shifted inward with respect to mass. 
Obviously this situation can only last because fresh material for 
the nucleosynthesis is transported downward by dredge-up. 

At the fifth TP (3 Mq) the core masses of the sequence 
with and without overshoot have not yet diverged too much 
due to the dredge-up difference between the two cases. How- 
ever, the differences due to overshoot are already well estab- 
lished. Not only is the mass-wise extent of the PDCZ larger 
with overshoot (Fig. [7|d) but in addition the convective insta- 
bility is longer (Fig.|lO|a). The He-flash convection zone devel- 
ops and disappears faster with overshoot. The overall duration 
of convective instability in the intershell is about 100 yr with 
overshoot and more than 150 yr for the case without overshoot. 

With overshoot the temperature at the bottom of the He- 
flash convection zone is larger while at the top of the con- 
vectively unstable zone the temperature is smaller (Fig.p^). 
Although the overall duration of the convective instability 
is larger without overshoot it is important to note that the 
high-T phase is much longer with overshoot. The duration 
of the high-T phase and the temperature reached during this 
phase is important for the analysis of the nucleosynthesis dur- 
ing the He-flash, in particular of the ^^Ne neutron source 
(^^Ne(a, n)^^Mg). Without overshoot the maximum temper- 
ature reached at the bottom of the He-flash convection zone 
during the fifth TP is 2.44 • 10^ K. With overshoot the tem- 
perature at the bottom of the He-flash convection zone exceeds 
2.5-10** K for 21.8 yr while T > 2.7- 10^ K for 9.3 yr. Over the 
following TPs the temperature in the He-flash convection will 
increase steadily. However, the general trend described above 
will remain preserved. Therefore, the exact mechanism of the 
s-process and also other aspects of the nucleosynthesis will be 
affected not only by the different intershell abundances due to 
overshoot but also because of different temperatures and time 
scales in the intershell during the He-flash. 



Fig. |10p displays the geometric evolution of the He-flash 
convection zone together with the position of the maximum 
energy generation by helium burning. The fact that the temper- 
ature at the top of the He-flash convection zone is lower with 
overshoot is closely related to the greater expansion of the in- 
tershell in this case. With overshoot the largest geometric extent 
of the He-flash convection zone is about 1.4 times the distance 
found without overshoot. 



5.2. The abundances in the intershell region 

Models with overshoot have larger intershell mass fractions of 
carbon and oxygen at the expense of hehum. In Fig. |l]the vari- 
ation of the intershell abundances from pulse to pulse is shown. 
These models with overshoot show qualitatively the same de- 
pendence on the pu l se number like models without ov ershoot 
( Schonberner, 197S] ; Boothroyd and Sackmann, 1988b ) but the 
quantitative abundances are very different. During the first TPs 
the ^^C and ^^O abundance increases strongly at the expense 
of ''He. After about six TPs the ^He abundance in the intershell 
reaches a minimum while '^^C and ^^O go through a maximum. 
After about five to ten additional TPs all abundances are level- 
ing at values which are similar to those found after the sec- 
ond or third TP. Quantitatively, overshoot with the efficiency 
/ = 0.016 increases the amount of carbon by about a factor of 
two and the amount of oxygen by a factor of 10 to 20. 

Test calculations in order to study the dependence of the 
intershell abundance on the efficiency of overshoot reveal that 
there is a relation between the overhoot parameter and the in- 
tershell abundance: larger overshoot leads to larger carbon and 
oxygen abundances and consequently smaller helium abun- 
dance. We have followed the evolution from the same start- 
ing model before the first TP over three TPs with different val- 
ues of /. In Fig. |l2| we display the abundance in the upper part 
of the intershell after the third TP. For / = the abundances 



are the same as in Schonberner (1979). With larger and larger 
overshoot efficiency, the amount of oxygen increases almost 
linearly. 

This correlation leads to a quite stringent constraint for 
the possible range for the overshoot efficiency / at the bot- 
tom of He-flash convection zone. We compare our intershell 
abundances with the surface abundances of the [WC]-type 
central stars of planetary nebulae (Koesterke and Hamann, 
1997; [Hamann, 1997|; pe Marco et al., 1998D and the PG 1 159 



stars (Dreizler and Heber, 1998; Werner et al., 1999). The 



spectroscopic abundance analysis shows that these stars are 
very carbon rich and also oxygen rich. Typically, one finds 
(He,C,O)=(0.50,0.33,0.17) (mass fractions) for PG 1159 stai's. 
Post- AGB stars become hydrogen-deficient because the inter- 
shell m aterial appears at the surface in the aftermath of a very 
late TP ( Herwig etal., 1999a ). Thus, the observed surface abun- 
dances of H-deficientpost-AGB stars must be interpreted as the 
intershell abundances of the progenitor AGB stars (Herwig and 
Blocker, 2000). Models without overshoot never show more 
than 2% of oxygen in the intershell and are thus unable to re- 
produce the large mass fraction of oxygen observed. This is a 
strong indication of the presence of at least some overshoot at 
the bottom of the He-flash convection zone. However, practi- 
cally no H-deficient post- AGB stars are known to show more 
than about 20% of oxygen^ 

The oxygen abundances shown in Fig.|l^ are an estimate 
of the surface oxygen abundance of a H-deficient post-AGB 
model evolved from the respective AGB model because the 
oxygen abundances after the third TP is already similar to the 



The exception, H1504+65, is free of H and He and has C ~ O. 
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later TPs (Fig. Any overshoot parameter much larger than 
/ = 0.03 would lead to H-deficient post- AGB models with 
too large an oxygen abundance. Therefore, observations - to- 
gether with the theoretical understanding of the evolutionary 
origin of H-deficient post-AGB stars - do constrain the over- 
shoot efficiency at the bottom of the He-flash convection zone 
to a narrow range of 0.01 < / < 0.03. 

5.3. The third dredge-up and its dependence on overshoot at 
the bottom of the He-flash convection zone 

One of the most important properties of TP- AGB models with 
overhoot is the efficient TDUP, which solves the problem of 
previous models to account for the typically low luminosities 
of carbon stars. In Table|l] and ^ the amount of dredged-up in- 
tershell mass and the related values for the dredge-up efficiency 
A = AMa/AA/hO are given for the 3 M© and 4 Mq sequence 
with overshoot, and these results will be further discussed in 
Sect.|. 

Here, we consider the dependence of A on the efficiency 
of overshoot. In particular it is important to clarify the differ- 
ent roles of overshoot below the bottom of the He-flash con- 
vection zone and below the bottom of the envelope convec- 
tion zone (see also Sect. 4-.1 ). For this purpose we have made a 
numerical experiment which involved the computation of one 
TP cycle with different values for the overshoot efficiency pa- 
rameter /. The computations were started immediately after 
the He-flash convection zone of the eighth TP had disappeared. 
Thus, the dredge-up episode following this TP has been com- 
puted with different / values at the envelope convection zone 
only, while during the preceding TP / was the same for all 
convective boundaries. The resulting values for A are shown 
in Fig.|l^ as open symbols and demonstrate the dependence 
of the third dredge-up on the overshoot efficiency below the 
convective envelope. Within the considered range of /-values 
practically no dependence of A on the overshoot efficiency ex- 
ists. Only without any overshoot (/ — 0) does dredge-up not 
occur for this TP. Even with / as small as 0.004 dredge-up is 
found with practically the same efficiency than with larger / 
(see Sect. \.\). Note that the value of A ~ 0.7 displayed for the 



open symbols in Fig. 13 is just the value found for the eighth 
TP of the 3 Mq overshoot sequence (Table[l|). If another TP 
had been chosen for this numerical experiment A would have 
had the correspondingly larger or smaller value. Fig.|l3| does 
not mean that models with overshoot applied only to the enve- 
lope convection always have A ~ 0.7. 

The test model sequences have then been evolved beyond 
the dredge-up episode after the ninth TP. The respective /- 
values have now been applied to the convective boundaries of 
the He-flash convection zone as well. The resulting A-values 
are represented by filled symbols in Fig. |l3| and show the de- 
pendence of the dredge-up on the overshoot efficiency at the 
bottom of the He-flash convection zone. The full symbol at 



/ = 0.016 shows A ~ 0.8, the value given in Table[l| for the 
ninth TP. / and A are correlated: larger overshoot leads to larger 
dredge-up. The relation defined by the full symbols is specific 
to the chosen TP of this sequence. However, the trend is gen- 
erally valid. A combination of effects is responsible for this 
correlation. With intershell overshoot the peak He-burning lu- 
minosities are substantially larger (Fig.^. At the 14'^ TP of 
the overshoot sequence (Mh — 0.658 Mq) the He-flash peak 
luminosity is Lho = 7.9 ■ IO'^Lq. At the 14*'' TP of the se- 
quence without overshoot (Afn = 0.692 Mq) the pea k lumi- 



nosity is Lhc = 2.4 • lO'^ Lq. Boothroyd & Sackmann ( |l988a| ) 



A = mass dredged-up after TP / growth of intershell by hydrogen 
burning during the preceding interpulse phase. 



found for a 3 Mq, Z=0.02 model sequence Lho — 2 • 10 Lq 
after about 20 TPs at a core mass A/h — 0.65 Mq. These 
models were computed without overshoot during the pre-AGB 
evolution and therefore have a lower core mass at the first TP 
A/hi ^ 0.55 Mq. 

Greater He-flash strength leads to a greater expansion and 
cooling of the upper layers of the intershell (Fig. |l^ and favors 
the occun-ence of the TDUP (Vrad oc T"*). In Fig.0 one 
can see that the modification of the intershell abundance during 
only one TP is considerable. The smaller He abundance with 
larger / leads to a larger opacity n which favors the dredge-up 
also because Vrad oc k. 

Thus the overshoot at the base of the He-flash convection 
zone strongly affects TDUP. In the previous section we have 
shown that / is constrained by the observational properties of 
H-deficient post-AGB stars. Therefore we conclude that the up- 
per limit of the efficiency of TDUP can be constrained from 
observational properties of post-AGB stars. 

6. Surface properties 

In the previous sections we have described the different mecha- 
nisms by which overshoot influences the model properties. We 
will now focus on the surface properties of the models with 
overshoot. 

Stellar parameters The comparison of the core mass- 
luminosity relation of models with and without e fficient 
dredge-up shows clear differences ( Herwig et al., 1998 ). While 
models without dredge-up follow a linear relation when the 
asymptotic regime has been reached, models with very efficient 
dredge-up (A ~ 1) continue to increase in luminosity even if 
the core mass no longer increases. Here, the continuing radius 
decrease leads - according to simple homology relations (Refs- 
dal and Weigert, 1970) - to an increase in luminosities. This ef- 
fect is most efficient over the earlier thermal pulses where the 
relative radius decrease per TP cycle is larger than after many 
thermal pulses when the core asymptotically resembles a white 
dwarf. The radius effect is responsible for the sub-luminous 
phase (compared to the luminosities expected from the core 
mass - luminosity relation) of the first few TPs, which is well 
known from any TP- AGB model sequence. However, the lumi- 
nosity evolution continues to be strongly coupled to the core 
radius evolution. This becomes apparent only if the core mass 
is prevented from growing continuously in accordance with the 
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Table 1. Dredge-up properties of the thermal pulses of the 3 M0 sequence with / = 0.016. The table gives the TP number, age 
(zero point set at first TP which shows dredge-up), peak luminosity of the helium burning shell Lho/ L©, amount of dredged-up 
hydrogen free material AMx/ Mq, core growth since last TP AA/h/ M©, dredge-up parameter A (see text), mass of hydrogen 
free core Afn, mass coordinate of top of pulse-driven convective zone Mpt, mass coordinate of bottom of pulse-driven convective 
zone Mpb and isotopic and elemental abundance ratios at the surface immediately after the TP. 



TP 


age 


i-Hc 




AA/h 


A 


Mh 


Afpt 






C/0 


16/17Q 


16/18Q 


yr 


10'^ Lq 




lQ--> 


Mq 


Mq 


Mq 


1 


-100079 


0.46 








0.63087 


0.62969 


0.60598 


20.4 


0.29 


194.7 


606.7 


2 


-50212 


1.17 




2.0 




0.63288 


0.63213 


0.60795 


// 


II 


// 


// 


3 


1 


2.13 


0.03 


3.1 




0.63593 


0.63547 


0.61003 


// 


II 


// 


II 


4 


48519 


2.63 


0.4 


3.7 


0.10 


0.63962 


0.63920 


0.61377 


20.5 


0.30 


194.8 


607.0 


5 


97563 


3.18 


1.0 


4.3 


0.23 


0.64352 


0.64315 


0.61819 


22.3 


0.32 


196.1 


610.9 


6 


146883 


4.63 


2.1 


4.7 


0.44 


0.64726 


0.64694 


0.62267 


26.9 


0.38 


199.8 


622.7 


7 


195842 


6.71 


2.7 


5.1 


0.53 


0.65030 


0.65002 


0.62688 


33.0 


0.45 


205.5 


640.4 


8 


243582 


8.80 


3.7 


5.3 


0.71 


0.65289 


0.65265 


0.63048 


41.4 


0.54 


213.9 


666.6 


9 


291084 


13.9 


4.6 


5.7 


0.82 


0.65482 


0.65461 


0.63328 


51.8 


0.64 


224.5 


699.7 


10 


338588 


19.0 


5.4 


6.0 


0.91 


0.65616 


0.65598 


0.63542 


63.8 


0.74 


236.6 


737.4 


11 


388433 


29.4 


6.1 


6.4 


0.98 


0.65717 


0.65702 


0.63701 


77.7 


0.86 


250.0 


779.1 


12 


438979 


42.9 


7.0 


6.8 


1.02 


0.65773 


0.65759 


0.63805 


93.0 


0.97 


263.7 


822.0 


13 


491151 


59.9 


7.5 


7.3 


1.04 


0.65804 


0.65804 


0.63871 


109.5 


1.08 


277.5 


864.9 



Table 2. Dredge-up properties of the thermal pulses of the 4Mq sequence with / = 0.016. For explanation of the displayed 
quantities see Table [l|. 



TP 


age 




AM;, 


AA/„ 


A 


Mh 


Mpt 


Mpt, 


12/13(-i 


C/O 


1G/17Q 


16/18Q 




lOti Lq 


10-^ 




M0 


Mq 


Mq 


1 


-80 


2.26 


2.39 






0.78259 


0.78239 


0.77172 


19.8 


0.31 


274.3 


602.7 


2 


11660 


3.70 


3.12 


1.91 


1.63 


0.78211 


0.78192 


0.77096 


22.4 


0.34 


275.1 


604.4 


3 


25990 


6.21 


3.94 


2.61 


1.51 


0.78160 


0.78144 


0.77017 


27.7 


0.42 


277.6 


610.0 


4 


42350 


11.5 


4.79 


3.17 


1.51 


0.78083 


0.78070 


0.76945 


34.9 


0.51 


282.7 


621.2 


5 


60290 


21.0 


5.41 


3.58 


1.51 


0.77962 


0.77953 


0.76865 


43.5 


0.62 


290.6 


638.6 


6 


79610 


38.1 


5.63 


3.93 


1.43 


0.77814 


0.77806 


0.76766 


52.8 


0.72 


300.4 


660.2 


7 


99500 


59.8 


5.90 


4.16 


1.42 


0.77667 


0.77661 


0.76656 


62.0 


0.82 


310.4 


682.1 


8 


120030 


92.6 


5.92 


4.45 


1.33 


0.77522 


0.77517 


0.76593 


71.2 


0.91 


320.5 


704.4 


9 


140790 


123.3 


0.61 


4.59 


1.33 


0.77389 


0.77384 


0.76472 


80.1 


0.99 


330.2 


725.8 


10 


161900 


142.6 


0.64 


4.74 


1.37 


0.77253 


0.77248 


0.76336 


89.0 


1.06 


340.5 


748.6 


11 


183480 


195.5 


0.65 


4.92 


1.32 


0.77102 


0.77128 


0.76187 


98.3 


1.13 


352.4 


774.6 



radius shrinkage, as in models with efficient dredge-up. In ad- 
dition, as Marigo et al. ( 1999 ) have pointed out, possibly up 
to on e third of the luminosity increase observed by Herwig 
et al. (1998) can be ascribed to the well known effect of the 



molecular weight increase in the envelope as a result of dredge- 
up of processed material. 

The rather short-term variation of the stellar parameters 
during and between the thermal pulses are affected as well. In 
general the surface luminosity reacts to the TP in the deep in- 
terior by a sudden ( A< ~ thermal timescale ~ 100 yr) and 
drastic luminosity decline of about 40% of the pre-TP lumi- 
nosity. It is followed (again on the thermal time scale) by an 
immediate luminosity jump which forms the more or less pro- 
nounced TP surface luminosity peaks (see Fig.^ at t = yr). 
The model sequence without overshoot shows only small TP- 
luminosity peaks after the first TPs with peak luminosities 



lower than the quiescent interpulse luminosity. Only after the 
third TP does the peak luminosity exceed the interpulse lumi- 
nosity and over the following TPs the ratio of peak and inter- 
pulse luminosity grows. The model sequence with overshoot 
shows very large peaks already after the first TP with a ratio of 
~ 1.6 between the peak and the interpulse luminosity which is 
gradually decreasing towards later TPs. At the fifth TP (Fig. |l^ 
the surface luminosity peak is similar for the two sequences. 
The actual local He-luminosity peak during a TP is followed by 
a much smaller secondary peak. As this secondary maximum 
decays in the intershell the surface luminosity does also de- 
cline to the second surface luminosity minimum at ^ 5000 yr 



in Fig. 15 Another difference caused by overshoot is a faster 
recovery of the H-burning shell after the TP (top panel Fig. p^] ). 
This is also reflected by the evolution of the surface parameters. 
For the model sequence with overshoot it takes only ~ 7600 yr 
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from the TP until a 70% level of the previous interpulse lumi- 
nosity at the surface has been regained while the sequence with- 
out overshoot requires ~ 11800 yr. This effect is even more 
pronounced at later TPs, however it is then intermingled with 
the effect of different core mass evolution due to the different 
efficiency of the dredge-up. 

Chemical evolution The chemical evolution of AGB models 
with overshoot is different compared to models without over- 
shoot due to four main effects: 

1 . The formation of aJ£C pocket in the upper part of the inter- 
shell region (Sect. 4.2). The consequences for the n-capture 
nucleosynthesis will be investigated elsewhere. 

2. Substantially larger dredge-up, which leads to more effi- 
cient transport of processed material to the surface. 

3. Higher temperature in He-flash convection zone during the 
TP(Fig.|l§. 

4. Depletion of helium and larger ^^O and ^^C abundance in 
the intershell due to deeper penetration of the bottom of the 
He-flash convection zone into the C/O-core. 

These effects are reflected by the chemical evolution at the stel- 
lar surface (Table |l] and ^. 



sumably also have experienced a different number of TPsJ^ 
However, apart from giants with C/0 « 1 and very low 
12q^13q x^iiQ (J-giants), a correlation between the two ra- 
tios is apparent. With respect to the model data this correla- 
tion must be interpreted as a time evolution. The model se- 
quences do generally follow the observed trend. However, the 
theoretical relation shows an offset to the observations and also 
a somewhat smaller slope. In addition the increment of the ra- 
tio from TP to TP is rather large and depending on mass loss 
several more TP are likely for these model sequences, lead- 
ing to very large ^^C/^^C ratios. A somewhat smaller over- 
shoot efficiency would result in a lower dredge-up efficiency 
and a smaller amount of dredged-up oxygen. While the first 
effect would decrease the steps at which the model sequence 
passes along the diagram, the latter would increase the slope of 
the theoretical relation somewhat. As another detail, the initial 
abundance ratios have a considerable influence on the theoreti- 
cal prediction. Most of the displayed stars are probably of lower 
mass and therefore their initial ^^C/^'^C ratio before the first TP 
was closer to ^ 10 (Gilroy, 1989) than around 20. This initial 
12C/13Q ratio determines the the slope of the ^'^GI^'^C - C/O 
relation as well, with a smaller initial ratio leading to a larger 
slope. Thus, models with initial masses around or below 2 Mq 
should show a larger slope due to this effect and should show a 
better agreement with the observations. 



C/O and carbon isotopic ratio Recumng and efficient dredge- Oxygen isotopic ratios Observationally the oxygen isotopic 

i3r ratios ^^0/"0 and ^^Or^O of AGB stars extend to much 
larger values than predicted by evolutionary models. Theo- 
retical models without overshoot cannot account for increas- 
ing ^^O/^^O and ^^0/^*0 ratios because dredge-up is only 
achieved with difficulty for the relevant masses. Even if this 
problem is circumvented by assuming dredge-up to be present 
then the dredged-up material leads only to a rather inefficient 
dilution of ^^O and ^^O in the envelope because the dredged- 
up material is depleted in these isotopes compared to the en- 
velope value. The ^^O abundance in the intershell of models 
without overshoot is similar to the envelope abundance and no 
modification can be expected there. 

1^0/170 and 



up transports primary C into the envelope after each TP. ■^C 
is not produced or destroyed in the stars of the mass range dis- 
cussed here. This leads to a steady increase of the ^^C/^^C ra- 
tio at the surface. At the same time also the C/O ratio increases 
from its initial value of ^ 0.3 and exceeds unity after a num- 
ber of TPs. Note that in calculations with overshoot ^^O is also 
mixed into the envelope at each dredge-up event. Because of 
this additional ^^O the C/O ratio grows more slowly with over- 
shoot than without overshoot if A where the same. However, 
this effect is more than compensated for by two other effects: 
the dredge-up itself is larger and the mass fraction of ^^C and 
^^O of the dredged-up intershell material is larger. This leads 
to a larger increase of the C/O ratio for each TP. For example, 
as can be seen from Table ^ the 4 Mq sequence becomes C-rich 
at the tenth TP contrary to the results of Forestini & Charbon- 
nel ( |1997| ) and Marigo et al. ( |1996| ) who both find that a 4 M© 
AGB star of solar metallicity does not become a C-star at all. 
Note that these authors assume dredge-up to operate and that 
HBB cannot prevent a 4 Mq star from getting C-rich. 

Fig. |6| compares observational data of the chemical com- 
position of AGB stars with the model data of the sequences 
with overshoot. The filled symbols in the lower left corner cor- 
respond to the first TPs. The observational data shows proper- 
ties of M- S- and C-giantQ stars with different masses (most of 
them smaller than those of the model sequences) which pre- 



1987) determined for 



Harris et al. ( 11985| ; 

16Q/18Q values between 600 and 3000. It is also interesting 
to note that the two ratios are correlated. Moreover, both ra- 
tios seem to be correlated with the neutron exposure associ- 
ated with the observed s-process abundances and also with the 
carbon abundance. Both, an increasing carbon abundance and 
an increasing neutron exposure are indicating more efficient or 
more frequent dredge-up because s-process elements and car- 
bon are not made in or at the bottom of the envelope. It is there- 
fore likely that the enormous oxygen isotopic ratios are related 
to the TDUP 



Forestini & Charbonnel ( |1997[ ) have assumed dredge-up 
and predict from an extrapolation of their models without over- 



M-giants 
giants:C/0 > 1.0 



C/O < 0.8, S-giants: 0.8 < C/O < 1.0, C- 
> 1.0. 



Giants without Tc have been excluded from Fig. |l^ because their 
abundance an omalv might b e due to transfer of mantle material from 
a companion (Habing, 1996). 
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shoot an increase of the ^^O/^^O ratio of 13% (2%) from ini- 
tially ^ 320 for the entire AGB evolution of their solar 3 M0 
(4Mo) case and almost the same values for the ^^0/^*0 ra- 
tio starting from initially ~ 600. The evolutionary effect in the 
overshoot model is significantly larger because the abundance 
of ^^O in the dredged-up material is five to ten times larger 
if overshoot is considered compared to models without over- 
shoot. Over the first dozen TPs available in Table ^ and ^ the 
oxygen isotopic ratios increase by about 44% (28%) for the 
3 Mq (4 Mq) case. Similar to what has been said previously a 
comparison of absolute numbers observed to the model predic- 
tions is misleading at this stage. The initial model ^^O/^^O ra- 
tios are probably too small, possibly due to a wrong ^^0(p, a) 
reaction rate. Also, most of the observed data will belong to 
less massive stars for which comparable model data do not yet 
exist. Clearly, this issue needs a more detailed investigation. 
However, it seems that overshoot can contribute to the solution 
of this problem. 



Magnesium isotopic ratios The evolution of the magnesium 
isotopic ratios ^"^Mg/^^Mg and ^^Mg/^^Mg are not only tracers 
of the TDUP but also of the temperature in the He-burning shell 
during the thermal pulses. While ^^Mg is not produced or de- 
structed in AGB stars the other Mg isotopes are produced dur- 
ing the TP by a captures of ^^Nc. The latter is abundant in the 
intershell as a result of two a captures on ^^N. Since the tem- 
perature at the bottom of the He-flash convection zone is larger 
if overshoot is present (Fig.jl^) the reactions ^^Ne(a, n)^^Mg 
and 22Ne(a, jf^Mg are more efficient. The 3 M© (4M0) se- 
quence accordingly shows a decrease of ^^Mg/^^Mg from 8.0 
(8.1) at the first TP to 4.9 (3.9) at the last TP computed. The 
24]V[g/26Mg ratio is 7.1 (7.0) for the 3 Mq (4 Mq) sequence at 
the first TP and decreases to 5.8 (5.3) at the last computed TP 
(see previous tables). Again, this decrease for the first dozen 
TP is stronger than that predicted by Forestini 8c Charbonnel 



( 1997) for the entire evolution. 



7. Conclusions 

We have studied how the evolution of AGB stars is affected by 
exponential diffusive overshoot at all convective boundaries. 
The main differences between models with overshoot com- 
pared to models without overshoot are 

1 . - a depletion of helium and enhancement of carbon and 

oxygen in the intershell abundance distribution which 
is caused by, 

- a deeper penetration of the He-flash convection zone 
into the C/0 core below (intershell or fourth dredge- 
up), 

2. - larger energy generation by He-burning during the ther- 

mal pulse and consequently, 

- higher temperature at the bottom of the He-flash con- 
vection zone and lower temperature at the top. 



3. - substantially increased dredge-up of material which has 

a modified abundance distribution compared to models 
without overshoot and as a result, 

- a different evolution of the surface abundances showing 
a stronger signature of nucleosynthesis and 

- a constant or decreasing core mass during some more 
advanced TPs, 

4. - a faster recovery of the hydrogen-burning shell after the 

TP and correlated with this, 

- a faster recovery of the surface abundance after the TP 
compared to models without overshoot, and 

5. the formation of a ^■^C and a ^^N pocket at the envelope- 
core interface at the end of the dredge-up. 

These properties of AGB evolution models improve the 
understanding of observationally obtained information about 
AGB stars. To date no obvious contradictions of observational 
properties and AGB models with overshoot have appeared. On 
the other hand, overshoot in AGB stellar models (in particular 
that at the base of the He-flash convection zone) has been an 
indispensable precondition for the starting point of H-deficient 
post-AGB model sequences which can explain the observed 
abundance ratios of these stars, in particular their oxygen abun- 
dance (Herwig et al., 1999a). Therefore, overshoot should be 
considered an important ingredient for the modeling of AGB 
stars. 
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Fig. 7. Evolution of the position and size in mass coordinates 
of several quantities at each TP as a function of time. Full sym- 
bols show the 3 Mq sequence computed with overshoot at all 
convective boundaries, open symbols represent the compari- 
son sequence without any overshoot and started with the same 
initial model before the first TP as the 3 Mq sequence with 
overshoot. Panel (a): H-free core at the TP, the increase by H- 
buming during the interpulse period is counterbalanced by the 
dredge-up after the TP for the overshoot sequence; (b): differ- 
ence of core mass at TP and smallest mass coordinate of con- 
vective envelope bottom achieved after the TP, negative values 
indicate DUP; (c): difference of core mass at TP and the largest 
mass coordinate of the PDCZ top achieved during the TP, (d): 
difference of top and bottom mass coordinate of PDCZ. 
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Fig. 8. Onset of He-flash convective instability in the inter- 
shell region at the first TP of 3 M0 sequence with overshoot. 
The two vertical lines show the upper (right) and lower (left) 
boundary of convection. 
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Fig. 9. The maximum helium-burning luminosities during the 
He-flash for the 3 Mq model sequence with (filled triangles) 
and without (open triangles) overshoot. Each symbol corre- 
sponds to one TP. The generally larger peak flash luminosities 
of the overshoot TPs are due to the extra-mixing below the bot- 
tom of the He-flash convection zone which leads to higher tem- 
peratures in the He-shell during the flash (Fig. 10). The relation 
runs vertically at M^i ~ 0.657 for the OV case because 
efficient dredge-up prevents further core mass growth (A ~ 1). 
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Fig. 10. Panel (a): The mass coordinates of the PDCZ bound- 
ary during the fifth TP of the 3.0 Mq with and without over- 
shoot. The time has been set to zero at the moment where the 
He-burning luminosity has reached the maximum. Panel (b): 
Temperature at the boundary of the He-flash convection zone 
of the same TP. Panel (c): Geometric position of the PDCZ 
boundary (CB) and of the maximum energy release by helium 
burning. 
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Fig. 11. Time evolution of abundances (mass fraction) in 
the upper part of the intershell shortly after each TP of the 
3 Mq and 4 M© model sequence with overshoot (compare with 
Boothroyd & Sackmann ( 1988a ), Fig. 9). Each mark indicates 
one thermal pulse starting at the first TP. The zero point for the 
3 Mq sequence has been set at the third TP where the TDUP 
stars. 
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Fig. 13. Variation of the dredge-up efficiency A with the over- 
shoot efficiency. The open symbols represent cases where the 
/ value has been changed just after the pulse-driven convective 
zone has disappeared after the eighth TP but before the bottom 
of the envelope convection has come close to the mass coordi- 
nate of the H-free core. Thus, in this case a different overshoot 
parameter is only effective at the bottom of the envelope dur- 
ing the dredge-up phase. The full symbols show the dredge-up 
efficiency after the following ninth TP which has been com- 
puted with the respective /-value also for the He-flash convec- 
tion zone. In this case a different overshoot parameter has been 
effective at the bottom of the He-flash convection zone. 
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Fig. 14. The intershell abundances after the ninth TP as a func- 
tion of /. While for Fig. |l^ the different overshoot values have 
been applied over all previous TPs here we have computed 
only the ninth TP with different overshoot values. The start- 
ing model has been a 3 Mq, / = 0.016 model before the ninth 
TP. This figure demonstrates how the abundances are modified 
due to overshoot during only one TP. 
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Fig. 12. Intershell abundances for different overshoot efficien- 
cies after the third TP. Test model sequences have been cal- 
culated from the same starting model before the first TP with 
different values of /. The abundances for the sequence with 
/ = are computed without overshoot. 
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Fig. 15. Luminosity of H-burning (top panel), effective temper- 
ature (middle panel) and surface luminosity (bottom panel) for 
one complete pulse cycle of the sequence with (solid line) and 
without (dashed line) overshoot respectively. Time has been set 
to zero at the fifth TP. 
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Fig. 16. Reproduction of i^c/isq versus ^"^C/^^C plot in 
Smith and Lambert (1990) (open symbols) together with the 
surface ratios of the 3 Mq (filled circles) and 4 Mq (fiUed pen- 
tagons) model sequence. See text for details. 



